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Abstract 

We describe two BQP-complete problems concerning properties of 
sparse graphs having a certain symmetry. The graphs are specified 
by efficiently computable functions which output the adjacent ver- 
tices for each vertex. Let i and j be two given vertices. The first 
problem consists in estimating the difference between the number of 
paths of length m from j to j and those which from i to j, where m 
is poly logarithmic in the number of vertices. The scale of the estima- 
tion accuracy is specified by some a priori known upper bound on the 
growth of these differences with increasing m. The problem remains 
BQP-hard for regular graphs with degree 4. 

The second problem is related to continuous-time classical random 
walks. The walk starts at some vertex j. The promise is that the dif- 
ference of the probabilities of being at j and at i, respectively, decays 
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with 0(exp(— fit)) for some [i > 0. The problem is to decide whether 
this difference is greater than aexp(— fj,T) or smaller than bexp(—fj,T) 
after some time instant T, where T is polylogarithmic and the differ- 
ence a — b is inverse polylogarithmic in the number of vertices. Since 
the probabilities differ only by an exponentially small amount, an ex- 
ponential number of trials would be necessary if one tried to answer 
this question by running the walk itself. 

A modification of this problem, asking whether there exists a pair 
of nodes for which the probability difference is at least aexp(— /iT), is 
QCMA-complete. 

1 Introduction 

Although it is commonly believed that quantum computers would enable us 
to solve several mathematical problems more efficiently than classical com- 
puters, it is difficult to characterize the class of problems for which this is 
expected to be the case. The class of problems which can be solved efficiently 
on a classical computer by a probabilistic algorithm is the complexity class 
BPP. The quantum counterpart of this complexity class is BQP, the class of 
problems that can be solved efficiently on a quantum computer with bounded 
error. Thus, the exact difference between the complexity classes BPP and 
BQP remains to be understood. 

An important way of understanding a complexity class is to find prob- 
lems which are complete for the latter. Meanwhile, some examples of BQP- 
complete problems are known [TJ |21 El HI E] • For the results presented here, 
the ideas of [5] are crucial. We rephrase the main idea and proof of this arti- 
cle since they provide the basis for the present work. The problem is stated 
in terms of sparse matrices. Following [H1I71IB], we call an N x N matrix A 
sparse if it has no more than s = poly log (N) non-zero entries in each row 
and there is an efficiently computable function which specifies for a given row 
the non-zero entries and their positions. 

We have formulated a decision problem concerning the estimation of the 
jth diagonal entry of the mth power of A and proved that it is BQP-complete. 
The proof that this problem is BQP-hard as well as the proof that it is in 
BQP both rely on the following observation. The jth diagonal entry of 
the mth power is the mth statistical moment of the probability distribution 
of outcomes when a measurement of the "observable" A is applied to the 
jth basis vector. Thus, the BQP-complete problem consists in estimating 
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the probabilities of the outcomes in an A-measurement. Moreover, we have 
shown in |Hj that the problem remains BQP-hard if it is restricted to matrices 
with entries ±1, 0. 

It is not surprising that questions related to decompositions of vectors 
with respect to eigenspaces of large matrices lead to hard computational 
problems and that quantum computer could be more powerful in dealing 
with such problems. After all, the dynamics of a quantum system is de- 
termined by the spectral resolution of its initial state. In this context, the 
interesting questions and challenges are to determine to what extent it is 
possible to identify natural problems related to spectral resolutions which 
are not mere reformulations of questions about the dynamical behavior of 
quantum systems. 

Calculating spectral decompositions with respect to large real symmetric 
matrices is certainly not a problem which only occurs in quantum theory. For 
example, methods in data analysis, machine learning, and signal processing 
rely strongly on the efficient solution of the principle component analysis [§] . 
However, data matrices typically occurring in these applications are in gen- 
eral not sparse and, even worse, their entries are usually not specified by an 
efficiently computable function. The data are rather given by observations. 

In our opinion, it is an interesting challenge to construct BQP-problems 
concerning spectra of sparse matrices having only and 1 as entries since 
problems of this kind are related to combinatorial problems in graph theory. 
It is quite natural to assume that a graph is sparse in the above sense; the 
efficiently computable function specifying the non-zero entries in each row is 
then a function which describes the set of neighbors of a given vertex and 
it is assumed that the number of neighbors is only polylogarithmic in the 
number of nodes. 

The paper is organized as follows. In Section |2] we present a quantum 
algorithm for measuring functions of observables. This algorithm uses the 
quantum phase estimation algorithm as a subroutine. We use this algorithm 
to show that the problems "Diagonal Entry Estimation" (Section |3J), "Differ- 
ence of Number of Paths" (Section^}, and "Decay of Probability Differences" 
(Section^ are all in BQP. Even though the first problem was already de- 
scribed in jS] we have first to extend these results in Section 01 since the proofs 
for the BQP-completeness of the other two problems are based on such an 
extension. In Section 0] we present a BQP-complete problem which (at first) 
strongly resembles the problem "Diagonal Entry Estimation". The impor- 
tant difference is that it deals with adjacency matrices (that is, 0-1-matrices) 
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instead of matrices with entries ±1,0. We describe how the problem is re- 
lated to mixing properties of a classical random walk with discrete time. In 
Section El we describe a BQP-complete problem concerning the mixing prop- 
erties of continuous-time classical random walks on regular sparse graphs. In 
Section El we describe a modification of the mixing problem which is QCMA- 
complete. 

2 Measuring continuous functions of observ- 
ables 

Before we explain what it means to measure functions of observables we 
have to introduce some terminology. The spectral measure induced by a 
self-adjoint N x iV-matrix A and a state vector 6 C N is defined as the 
measure on R supported by the set of eigenvalues Xj of A with probabilities 
Pj := (ip\Qj\ip), where Qj is the spectral projection corresponding to Xj. 

Let / : R — > R be a function. The spectral measure induced by f(A) and 
\if}) is then the measure on R supported by the values f(Xj) with probabilities 
Pj, where pj and Xj are as above. 

We define measuring a function / of an observable B in the state as 
a quantum process which allows us to sample from a probability distribution 
that coincides with the spectral measure induced by f(B) and Note that 
when / is the identity function, then this corresponds what is considered to 
be a von-Neumann measurement of an observable in quantum mechanics. 
As in [J] the main tool for implementing such measurements is the quantum 
phase estimation algorithm ^U|. The idea is that for each observable B with 
< 7i we can implement the measurement of B by applying quantum 
phase estimation to the unitary exp(lB), i.e., the map that describes the 
corresponding dynamics according to the Hamiltonian —B. For the class of 
self-adjoint operators considered here an efficient simulation of the dynamics 
is indeed possible since they are sparse [HI E] • By applying the procedure 
to a given state vector \ip) we can sample from the spectral measure induced 
by B and 

The main statement of this section is that we can also sample from the 
spectral measure induced by f{B) since we can implement such von-Neumann 
measurements of B and then apply / to the outcomes. Hence, no implemen- 
tation of exp(if(B)) is needed (this is important because, for sparse B, the 
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matrix f{B) will in general not be sparse and therefore it is not clear how 
to simulate the corresponding dynamics). In particular, we can estimate the 
expectation value of a given / with respect to the spectral measure. How- 
ever, for our proofs it is essential how the accuracy behaves when certain 
smoothness assumptions are posed on /. This is made precise in Lemma [0 
But first we have to recall how quantum phase estimation works jTUJ . We 
add a p-qubit ancilla register to the n qubits on which V := exp(iB) acts. 
Then we replace the circuit implementing V with a controlled-V gate by 
replacing each elementary gate with its controlled analogue (note that this 
cannot be done in a black-box setting). The ancilla register is initialized to 
the equally weighted positive superposition of all binary words, then we use 
copies of the controlled gate to apply the 2 J th power of V to the n qubits 
if the jth ancilla qubit is in the state |1). Finally, we apply the inverse 
Fourier transform to the ancilla register and measure the ancillas in the 
computational basis. Given that the n-qubit register was in an eigenstate of 
V with eigenvalue exp(i2ir(p) the obtained binary word a G {0, 1, . . . 2 P — 1} 
provides a good estimation for 27i(p in the following sense: 

Pr(\<p - a/2 p \ < ij) > 1 -6, (1) 

where 9, rj > and the number of ancillas is given by jTTj 

p:= riog(l/77)l + riog(2 + (l/(20))l. (2) 

Here, the distance \<p — a/2 p \ is to be understood with respect to the cyclic 
topology of the unit circle which identifies if = 1 with ip — 0. In contrast to 
the convention in JT] we reinterpret throughout the paper a as phases x in 
the interval [— ir, tt) by defining x := a2n/2 p if a < 2 P ~ 1 and x := a2ir/2 p — 27i 
otherwise. We refer to x as the outcome of the measurement procedure. Then 
we have: 

Lemma 1 (Continuous Functional Calculus) 

Given a self- adjoint n-qubit operator B with < 1 and a quantum circuit 
U with || {7 — exp(iB)|| < 5 for some 5 > such that the decomposition of U 
into elementary gates is known. Let be an n-qubit state whose decompo- 
sition into B- eigenvectors contains only eigenvectors whose eigenvalues are 
in the closed (but possibly infinite) interval I . Let f : I — > R be a Lipschitz 
continuous function with constant K, i.e., \f(x) — f(y)\ < \x — y\K for all 
x,y e I. 
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First, we apply the phase estimation (with unitary U) to , where the 
number p of ancilla qubits is given by eq. HJ). Next, we apply f to the outcome 
x if x & I, otherwise we replace x by the value in I closest to it and apply 
f to the latter. The expectation value of the random variable f(X) in this 
experiment, denoted by E^\(f(X)), satisfies 

\E m (f(X)) - W(B)\iP)\ < (26 + 52 p+1 ) H/IU + 2-kK V , 

where r/,9 and p are related to each other as in eq. (OJ) and \\f\\oo denotes the 
norm given by the supremum of \f(x)\. 

Proof: We first assume we could implement V := exp(iB) instead of its 
approximation U. In this case, if we apply phase estimation to an eigenvector 
\ipj) of B with eigenvalue Xj, then the outcome x is likely to be close to A-,- 
in the sense that 

PidXj -x\< 2ttt]) >l-6. 

This follows from ineq. and the fact that the values a/2 p have to be 
rescaled by 2tt. We conclude 

\E^ } (f(X)) - f(Xj)\ < 20H/IU + 2nK V . (3) 

The second term on the right of ineq. (jSJ) corresponds to the case that x 
deviates from Xj by at most 27rrj (note that the probability for this event 
is not decreased by replacing all outcomes lying outside from / with the 
closest values in /). Then the Lipschitz condition ensures that the error of 
f(x) is small. If the error of x is large (which happens with probability at 
most 9), the error of f(x) is still upper bounded by 2 ||/||oo- If we replace 
the eigenvectors \ipj) with some vector = Cj\i[)j) , the distribution of 
outcomes in the phase estimation procedure is a mixture of the distributions 
for each with weights \cj\ 2 . This is easily checked by analyzing the 
standard phase estimation. Thus we obtain 

\E m (f(X)) - < 20H/IL + 27rKr]. (4) 

from ineq. (J3J) by convexity arguments. 

Now we take into account that we have only an approximation U of 
V. Since the procedure uses 2 P+1 — 1 copies of the controlled-t/ respective, 
controlled-^, the norm distance between the phase estimation circuits is less 
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than 2 P+1 5. Let q(a) and q(a) denote the probability of obtaining the result a 
when the phase estimation is implemented with U and V, respectively. Then 
the / 1 -distance between the measures q and q satisfies (for technical details 
see jS]) 

\\q-q\\i:=J2M a )-^\^ 2P+2S - 

a 

The error in the expectation value of f(X) caused by applying / to outcomes 
x computed from a when a is sampled according to q instead of q is hence 
less than 2?+ 2 5 WfW^. 

Putting everything together we obtain the desired bound 

\E m {f(X)) - (i/>\f{B)\i/>)\ < (26 + 2"+ 2 6) ||/|U + 2tt^. (5) 

Since we have shown that the expectation value of the values f(x) in the 
procedure in Lemma Q is close to the expectation value (ip\f(B)\tfj) we can 
use the average over the values f(x) after a small number of runs as a good 
estimate for the desired expectation value. The following lemma states that 
we obtain an efficient procedure for giving an estimation up to any desired 
accuracy which is inverse polynomial in n. 

Lemma 2 (Estimation of Expectation Values) 

Given an efficient algorithm to simulate exp(iB) on n qubits with \\B\\ < 1 
in the sense that the resources to obtain a circuit U with \\U — exp(iB)|| < 5 
scale polynomially in n and 1/5. Given a state whose decomposition 
into B- eigenvectors contains only those with eigenvalues in the interval I. 
Let f be as in Lemma[I\ Then we can estimate (ip\f(B)\ip) up to an error 
£ ■ (||/||oo + K) with probability at least 1 — a such that the time and space 
resources are bounded by a polynomial in n, 1/e, and log(l/a). 

Proof: The empirical average of the values f(x) obtained after a few runs 
converges exponentially fast to the expectation value E\^(f(X)). Using Ho- 
effding's bound [12] and the fact that f(X) is a random variable with bounded 
range (for details cp. also one can easily show that the required number 
of runs for estimating then expectation value up to an accuracy e scales in- 
verse polynomial in e and polylogarithmically in I /a when this error should 
be guaranteed with probability 1 — a. 

It remains to bound the resources required to ensure that 

\E m (f(X))- W(BM)\ < e- (\\f\\oo + K). 
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Taking into account that the number of required ancilla qubits p is only 
logarithmic in 1/9 and 1/rj (see eq. (0)), ineq. (j3J) shows that the error is 
bounded by (||/||oo + K) times a polynomial in 1/9, 1/rj, 1/5 and n. 

In order to make the overall error smaller than e ■ (\f(I)\ + K) we ensure 
that each of the three terms 29 ||/||oo 5 2 P+1 6 ||/||oc 5 an d 2nKr] on the rhs. of 
ineq. © is at most e • (1/(7)1 + K)/3. 

To this end, we choose rj = e/(67r) and 9 := e/6. Putting these values 
into eq. (0) we obtain the number of required ancilla qubits p. We choose 
5 such that S2 P+2 < e/3. Note that 2 P is polynomial in 1/e, so that 1/5 is 
polynomial in 1/e and n. The number of controlled-f/ circuits is a multiple 
of 2 P and, thus, only polynomial in 1/e. 

Lemma|21shows that there is an efficient quantum algorithm for measuring 
functions of observables provided that the demanded accuracy is not too high. 
In the following sections we introduce several problems and prove that they 
are BQP-complete. The proofs that they are in BQP rely on the quantum 
algorithm for measuring functions of observables. 

3 Diagonal Entry Estimation 

In this section we extend the results of [3] saying that the estimation of di- 
agonal entries of A m for a sparse matrix with entries ±1, is BQP-complete. 
Here we argue that the construction in j3] shows also that the problem re- 
mains BQP-hard if when further restricting to matrices with only 4 non-zero 
entries. This is made precise in Lemma El 

Furthermore, Section El makes it possible to provide a unifying picture of 
the results in [S] and our results presented in Sections 0] and El since for all 
three problems we use spectral measures of functions of observables to show 
that they are in BQP. 

We first state the formal definition of the decision problem related to the 
estimation of the diagonal entries of powers of sparse matrices. 

Definition 1 (Diagonal Entry Estimation) 

Given a sparse symmetric N x N -matrix A with real entries, an integer 
j G {1, . . . , iV} ; and a positive integer m = polylog(iV) ; estimate the diagonal 
entry (A m )jj in the following sense: 
Decide if either 

{A m ) 3j >g + tb m 
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or 

(A m ) J3 <g-eb m , 

for given g G [— b m , b m ] and e = l/polylog( N) , where b is an a priori known 
upper bound on the operator norm of A. 

We showed in 

Theorem 1 The problem "Diagonal Entry Estimation" is BQP-complete. 
This remains true if the matrices A have only ±1,0 as entries. 

We recall the formal definition of the complexity class BQP |13| . 

Definition 2 (The class BQP) 

A language L is in BQP if and only if there is a uniformly generated family 
of quantum circuits Y r acting on poly(r) qubits that decide if a string x of 
length r is contained in L in the following sense: 

F r |x,0) =a x , |0)® |^x,o) +a x ,i|l) ® |^x,i) (6) 

such that 

1- |a x ,i| 2 > 2/3 if-x e L and 

2- |a Xi i| 2 < 1/3 ifxg_L. 

Equation (QJ) has to be read as follows. The input string x determines the 
first r bits. Furthermore, I additional ancilla bits are initialized to 0. After 
Y r has been applied we interpret the first qubit as the relevant output and 
the remaining r + I — 1 output values are irrelevant. The size of the ancilla 
register is polynomial in r. 

To show that "Diagonal Entry Estimation" is in BQP we simply apply 
the quantum algorithm for measuring functions of the observable A/b to the 
state \j), and the function / : [— 1, 1] — > R with f(x) := x m . Then we can use 
Lemma El to show that we can efficiently estimate (j\A m \j) with sufficient 
accuracy. To see this, we observe that the Lipschitz-constant K satisfies 
K < m and we have ||/||oo — 1- Hence we can efficiently estimate the 
diagonal entries of (A/b) m up to e (m + 1) for any desired inverse polynomial 
e. With e' := e (m + 1) we have an accuracy e'. Thus, we may achieve 
an accuracy of e' b m for the diagonal entries of A m for any desired inverse 
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Figure 1: Circuit £/ constructed from the original circuit Y r . Whenever the answer 
of the BQP problem is no, the output state of U is close to the input state |x, 0) = 
\x±, . . . , x n , 0, . . . , 0). Otherwise, the state |x, 0) is only restored after applying U 
twice. 

polynomial e'. By Definition ^ this is sufficient to solve "Diagonal Entry 
Estimation" . 

The idea of our proof of the BQP-hardness is based on an encoding of the 
quantum circuit which solves a given BQP-problem into a self-adjoint opera- 
tor A such that the spectral measure induced by A provides the information 
on the solution. To this end, we assume that we are given a quantum circuit 
Y r which decides whether a string x is in the given language L in the sense 
of Definition EJ 

In some analogy to ^2J|3] we construct a circuit U which is obtained from 
Y r as follows: first apply the circuit Y r , then apply the o^-gate on the output 
qubit, and finally apply the circuit Yj . The resulting circuit U is shown in 
Fig. We denote the dimension of the Hilbert space U acts on by N. 

Let U be generated by a concatenation of the M elementary gates 

Uq, . . . , Um-i ■ 

We assume furthermore that M is odd, which is automatically satisfied if we 
decompose Yj in analogy to Y r and implement a cr^-gate between Y r and Yj . 
We define the unitary 

Af-l 

W:= |J + l)<*|®E/i, (7) 

1=0 

acting on C M (g> C^. Here the + sign in the index has always to be read 
modulo M. 

Now we define the self-adjoint operator 

A:=^(W + W ] ). (8) 
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Then we apply an A-measurement to the state 



|«x) := |0) ® |x,0). 

Since |s x ) is a basis state, we can choose j such that \j) = |s x ). To understand 
why |a | influences the measurement statistics it is useful to consider the 
extreme cases = 0, 1. One may check that the repeated application of 
W to \j) leads for a = to an orbit (of mutually orthogonal states) that is 
periodic after M steps. The spectral measure induced by \j) and W is hence 
the uniform distribution over the Mth roots of unity. For |cni| = 1 the orbit 
is M-periodic up to the phase —1. This leads to a uniform distribution over 
the roots of unity reflected at the imaginary axis. As shown in jH] in detail, 
the spectral measure induced by \ j) and A depends therefore on a > a i m the 
following way: 

The spectral measure is the convex sum 

P:= \a \ 2 pW + \ ai \ 2 P^, (9) 

where P^ is supported by the values \\ := cos(27r//M) with probabilities 
p/ 0) = 1/M for I = and P z {0) = 2/M for / = 1, . . . , M - 1/2. The measure 
is obtained from P^ by a reflection at the origin and is hence supported 
by the values := cos(vr(2/ + 1)/M) with probabilities -Pffl-ij/a = 1/M 

and Pi = 2/M for I < (M — l)/2. The measure P can also be written as 

P(A) := ^(|«o|\ A (o)P/ 0) + laxW^Ph ■ (10) 
i 

Then we obtain the jth diagonal entry of A m as the mth moment of the 
measure P: 

= a-N')j:(Af)"if' + w'j:(A?')"if>. (id 

Z I 

Since the possible range of \a\\ is determined by the solution of the decision 
problem the latter also determines the range of the possible values for (A m )jj. 
One checks easily that (A m )jj changes sufficiently when changing |ai|, i..e, 
an estimation of (A m )jj makes it possible to decide whether xGL. 
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To see that the proof of BQP-hardness works also with matrices having 
only ±1, as entries we make the following observation. Hadamard-gates and 
Toff oli- gates form a universal set. In [S] we have instead chosen a universal set 
containing only Hadamard gates and gates which are concatenations of one 
Hadamard gate with a Toffoli gate. Then the above unitaries Ui have only 
±1/a/2 as entries. Since the decision problem Diagonal Entry Estimation is 
formulated in a scale invariant way, we may rescale A by y/2 and obtain only 
entries ±1, 0. 

In the following we will need the following extension of the results in [3] : 

Lemma 3 The problem Diagonal Entry Estimation remains BQP-hard if A 
contains only entries ±1,0 and the number of non-zero entries in each row 
of A is 4. 

The proof is obvious after observing that A as constructed above has 
exactly 4 nonzero entries in each column and, similarly, in each row. Each 
gate Ui has exactly two non-zero entries since it is the embedding of either 
a Hadamard gate or a Hadamard gate combined with a Toffoli gate. One 
checks easily that W has also two non-zero entries in each row leading to 
four entries in W + . 

4 Difference of Numbers of Paths 

In this section we show that a modified version of Diagonal Entry Estimation 
is still BQP-complete when restricting to matrix with entries 1, 0. This makes 
the problem more combinatorial than the original one. In other words, we 
introduce a problem concerning regular sparse graphs and show that it is 
BQP-complete. A graph G is called regular if the degrees of all its vertices 
are equal. We call G sparse if its adjacency matrix A is sparse, i.e., for each 
vertex the number of adjacent vertices is polylogarithmic in the number of 
nodes and there is an efficiently computable function that computes the set 
of adjacent vertices for each vertex. We define the following problem: 

Definition 3 (Difference of Numbers of Paths) 

Let A be the adjacency matrix of a regular sparse graph G on N vertices and 
m a positive integer m = polylog(iV). Let q and r be two vertices such that 
there exists an automorphism of G which exchanges q and r. Set 

A£0 :=(A m ) qq -(A m ) qr . 
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Decide if either 

A^>g + eb m 

or 

A^<g-eb m , 

for given g e [— b m , b m ] and e = l/polylog(iV) ; where b is an a priori known 
upper bound on 

sup|A^| 1/n - 

n 

Theorem 2 The problem "Difference of Numbers of Paths" is BQP-complete. 

We prove this theorem by showing that this problem is BQP-hard and is 
contained in BQP. The proof of BQP-hardness is based on a reduction to 
the problem "Diagonal Entry Estimation". An efficient quantum algorithm 
for estimating the difference of number of paths is obtained by applying the 
algorithm for measuring functions of observables. 

4.1 "Difference of Numbers of Paths" is BQP-hard 

We show that "Diagonal Entry Estimation" reduces to "Difference of Number 
of Paths". Let A be an N x N matrix with entries ±1, as in Definition ^ 
The reduction relies on the conversion of A into a 0-1-matrix A acting on a 
space of dimension D := 2N such that the restriction of A to some invariant 
iv-dimensional subspace is unitarily equivalent to A. 

The conversion works as follows. Each — 1 in A is replaced by the Pauli- 
matrix a x and each 1 by the identity matrix 1 2 = o~ 2 x . Each is replaced by 
the zero matrix O2. More formally, the ±l-0-matrix 

AT-l 
i,j=0 

is converted to the 0-1-matrix 

N-l 

A:= £ Ty ® 

i,j=0 

where 7^ := a x if = —1, 7^ := 1 2 if a^- = 1, and 7^ := 2 if = 0. 
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We now prove that the restriction of A to some iV-dimensional subspace 
is unitarily equivalent to A. To this end, we start by introducing the vectors 

1^) :=^(|0)±|1))GC 2 . 

On the subspace spanned by \<p~) the matrices 7^ act by multiplying the 
vectors by the scalars a^, i.e., 

7«I0~) = a ij\^) ■ 

Similarly, on the subspace spanned by they act by multiplying the vec- 
tors by the scalars |ay| = a 2 -, i.e., 

Consequently, we have the direct sum decomposition 

A = |<r>(<TI ® A + |0 + >(0 + | <gi(A*A), 

where A* A is the Hadamard product (or entry- wise product) of A with itself 
and has therefore only 0, 1 as entries. In the following, it will be important 
that A is invariant under conjugation with a x ® 1^. This invariance follows 
directly from the direct sum decomposition and the fact the vectors are 
eigenvectors of a x with corresponding eigenvalues ±1. 

Assume we want to estimate (A m )jj. Let \j) be the jth basis vector in 
C N . Then we define vectors in C 2 ® C N = C 2N by \ip~) := \<p~) ® 
\q) '■= |0) ® |r) := |1) ® and label the corresponding vertices by q 
and r, respectively. We have then = (\q) + \r))/y/2. Note that the 
automorphism a x ® 1 exchanges the vertices g and r. We reproduce diagonal 
entries of A m from both diagonal and off-diagonal entries of A m as follows: 

(A m )n = (^-\A m \^-) 
= V2(ii-\A m \q) 

= {A m ) qq - {A m ) qr = . (12) 

The second equality follows from the fact that the space \<fi~) ® C N is A- 
invariant. Thus, only the component of \q) lying in this subspace is relevant. 

The terms in eq. (112)1 correspond to the number of paths of length m 
from q to q and from g to r, respectively. In other words, the problem of 
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estimating diagonal entries of A m reduces to estimating the number of paths 
of length m from q to q and q to r in the graph defined by A, where q and r 
have to be chosen suitably. 

Taking into account that also off-diagonal entries of powers of sparse 
matrices can be estimated efficiently 0, it seems as if we had already shown 
that estimating entries of powers of adjacency matrices is BQP-complete. 
However, this is not the case. There is a subtle but very important issue 
concerning the precision of this estimation. The accuracy in Definition ^ is 
given by eb m , where b is an a priori known bound on the operator norm of 
A. In general, this bound will not be a bound on the norm of A. Using 
the direct sum decomposition of A, we see that the norm of A is equal to 
max{||A||, and is smaller or equal to max{6, b 2 } as < 

[IB] . On the one hand, if we estimate the entries of A m up to an accuracy 
specified in terms of the norm of A, then these values are not good enough to 
obtain a sufficiently precise estimate of the diagonal entry (A m )jj. Assume 
for instance, that c||A|| = ||A|| for some constant c > 1. Then the entries 
of A can only be estimated with an accuracy that is an inverse polynomial 
multiple of c m ||y4|| m , which is exponentially worse than the required accuracy 
e||A|| m . On the other hand, the operator norm \\A\\ of A is not a natural 
quantity to define the required accuracy of the estimation of the number of 
paths in A. 

For these reasons the estimation accuracy of the decision problem "Dif- 
ference of Numbers of Paths" is not (and should not be) formulated in terms 
of the norm of A. Instead, it is specified in terms of a quantity which is 
sufficiently small, making it possible to reduce Diagonal Entry Estimation 
to the former. For doing so, we recall that for each nGN the nth diagonal 
entry of A n can be expressed in terms of due to 

(A n h = A<? . 

We have therefore 

sup = sup |(^-|i n |VO| = sup|(A n y < \\A\\ n . 

neN neN neN 

This proves that if we can solve "Difference of Number of Paths" then we 
can also decide between the two cases in Definition^ with b := \\A\\. 

It remains to show that A is the adjacency matrix of a regular graph. 
We first recall Lemma El By replacing the values ±1 with the 2x2 identity 
matrix and the matrix a x , respectively, we obtain a matrix with four entries 
in each row, as well. Hence all vertices of the graph have degree 4. 
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4.2 "Difference of Number of Paths" is in BQP 



To solve the problem "Difference of Numbers of Paths" on a quantum com- 
puter we observe that ApP is an expectation value of a random variable 
which occurs in an appropriate measurement procedure. In analogy to the 
preceding paragraph, we define 

\r) :=^(|?> + |r», 

with the essential difference that the nodes q and r are a priori given and 
not derived from setting \q) := |0) (g) \ j) and |r) := |1) (g> \j). Then we have 
nevertheless 

A™ = (^\A m \^-), (13) 

which follows here from the symmetry defined by the given graph automor- 
phism. Without this symmetry, we had 

{r\A m \r) = i(A£> - 2Aj?> + A&>) , 

but here we have A^ = Affl. Due to the promise on the growth of A$ 
the decomposition of in A-eigenvectors contains only those eigenvectors 
whose corresponding eigenvalues satisfy Xj G [—b, b] . We may rescale A 
by some factor d such that its norm is at most 1. On the rescaled interval 
[—b/d, b/d] we define the function f(x) := x m having Lipschitz-constant K < 
m(b/d) m and sup-norm ||/||oo = (b/d) m . Then the achievable estimation error 
for (^~\(A/ d) m \il)~) = Aq-f* jd m is an inverse polynomial multiple of (b/d) m , 
i.e., we can estimate A^ up to any desired inverse polynomial multiple of 
b rn . 



4.3 Interpretation of difference of number of paths in 
terms of discrete-time random walks 

Note that there is also an alternative interpretation for the difference of the 
number of paths. Since A describes a regular graph of degree 4, we obtain a 
doubly stochastic matrix by 
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It describes a classical random walk (in discrete time) on the corresponding 
graph. Then the entry {A m )ij is the probability of reaching position j after 
m steps given that the initial position was %. For increasing m, the differ- 
ence (A m )ii — (A m )ij is therefore directly linked to the decay of probability 
differences, i.e., to mixing properties of the random walk. 

5 Decay of Probability Differences 

We can also formulate a BQP-complete problem in term of continuous-time 
random walks on regular sparse graphs. The dynamics of a random walk is 
generated by the Laplacian of the underlying graph. Consequently, the time 
evolution of an initial probability distribution is determined by its spectral 
resolution into the eigenvectors of the Laplacian. 

For regular graphs, the Laplacian coincides with the adjacency matrix up 
to the negative sign and an additive multiple of the identity matrix. Hence, 
the spectral resolution with respect to the adjacency matrix (being crucial in 
the preceding section) also determines the behavior of the continuous random 
walk. We introduce some notation to make these statements more precise. 

Let A describe a regular graph on iV vertices with degree d. We define 
its Laplacian by jTo] 

C:=dl-A. 

Let p(t) = (pi(t), . . . ,PN(t)) be a probability distribution where Pj(t) for 
j = 1, . . . , N is the probability of being at vertex j for j = 1, . . . , N. The 
Laplacian L defines a continuous-time (classical) random walk by 

p{t) := e- ct p, 

where p = p(0) is the initial probability distribution. 

Our BQP-complete problem considers the following question. Assume 
that the random walk starts at vertex q. Given a second vertex r, determine 
the decay of the difference between the probabilities p q (t) and p r {t). 

For infinite time, only the smallest eigenvalue whose eigenvector is con- 
tained in the initial state p determines the exponent dominating the decay. 
For our specific problem (where only the difference between the probabilities 
of two vertices is considered), only the spectral decomposition of |g) — |r) is 
relevant. For a random walk which started at vertex q, the difference between 
p q (t) and p r (t) is given by 

c, r (i):=(e-%-(e-%, (14) 
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where we use {W) qr to denote the entry of an arbitrary matrix W in row q 
and column r. 

Our BQP-complete problem is defined as follows. 

Definition 4 (Decay of Probability Differences) 

Let c qr (t) be defined as in eq. where £ is the Laplacian of a regular 

graph G on N vertices with degree d G 0(polylog(iV)) such that there is an 
efficiently computable function that determines the set of adjacent vertices 
for every vertex. Let fi > and b < a < 1 be two positive numbers with 
1/a, 1/6, l/(a — b) G polylog(iV) and T G polylog(iV) be some time instant. 
Let q, r be two vertices for which there is an automorphism of G exchanging 
q and r. Given the promise that 

c qT (t) G Ofe"" 4 ) , 

decide if either 

C qr (T) > fle"^ 

or 

c qr (T) < be~^ T . 

Note that in both cases, the exponent of decay may be the same but only 
the constants a and b differ. 

Theorem 3 The problem "Decay of Probability Differences" is BQP-complete. 

We prove this theorem by showing that the problem is BQP-hard and is 
contained in BQP. 

5.1 "Decay of Probability Differences" is BQP-hard 

To see that every problem in BQP can be reduced to "Decay of Probability 
Differences" we consider the adjacency matrix A obtained from the matrix 
A as described in Section 0] Adopting the notation used there, we have 
\q) := |0) ® \j) and |r) := |1) ® \j) and label the corresponding vertices 
by q and r, respectively and obtain in straightforward analogy the following 
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equality: 

c qr (t) = (q\e- Ct \q) - (q\e- ct \r) 

= V2(ij-\e- ct \q) 

= (ij-\e- ct \ij-) 

= (^-\e- {4 - l2N - A)t \^-) 

= (j\e-^»- A »\j) . (15) 

We need some properties of the spectral measure induced by A and \j) to 
prove that "Decay of Probability Differences" is BQP-hard. They are ob- 
tained from Section |3] after taking into account the rescaling factor a/2 (see 
the end of the section) and read: 

1. the largest eigenvalue in the support of the spectral measure is \/2 and 
the second largest is \/2 cos(7r/M), and 

2. the probability of obtaining the eigenvalue \/2 is equal to \oto\ 2 /M. 

Set /i := 4— a/2 and v := 4— a/2 cos(7r/M) which are equal to the smallest and 
second smallest eigenvalue, respectively, that occur in the spectral measure 
induced by 4 ■ In — A and Using these eigenvalues, eqs. (|T3|l imply 

c qr (t) > i^e"" 4 

and 



c qr {t) < ^e~^ + e- vt 



| «o| 
M 

Choose the time instant T := ln(6M)/(z/ — /i). In this case, the bounds are 

cAT) > ^e-^ 

and 

Note that T increases only polynomially with M since v — fi scales inverse 
polynomially with M. Taking into account that |ao| 2 <l/3ifxGL and 
l^o 1 2 — 2/3 otherwise, we have for x£l 

c,AT) < ^ , 
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and for x L 



CqriT) > 




This proves that we can reduce the question whether x e L to the estimation 



5.2 "Decay of Probability Differences" is in BQP 

In strong analogy to Section the quantity of interest, which is here c qr (T), 
can be written as an expectation value 



where the last equation follows (in analogy to Subsection 14.2)1 from the sym- 
metry of the graph and not, as in the preceding paragraph, from setting 
\q) := |0) <S> The last expression in eq. (|TT))l is an expectation value of a 
random experiment where the observable C is measured and the result A is 
converted to exp(— Xt). To show that "Decay of Probability Differences" is in 
BQP we rescale C to C/b where b is a polynomial upper bound of the norm 
of C. We have, for instance, ||£|| < 2d since the norm of A is given by its 
degree d . From the promise saying that the probability differences decay 
with 0(exp(— fri)) it follows that the decomposition of \ip~) contains only C- 
eigenvectors with eigenvalues in [/x, oo], i.e., we have only eigenvalues of C/b 
in [fi/b, oo]. We define / : oo] — ► R + with f(x) := exp(— xbT). Its Lip- 
schitz constant is K = 6Texp(— fiT). Moreover, we have ||/||oo — exp(— fiT). 
Lemma El guarantees that we can efficiently estimate the expectation value of 
f(C/b) up to any desired inverse polynomial multiple of (bT + 1) exp(— fiT). 
Since b and T are both polynomial, we can obtain an accuracy being any 
desired inverse polynomial multiple of exp(— fiT). This is sufficient to solve 
"Decay of Probability Differences" . 

6 A QCMA-complete mixing problem 

The complexity class QCMA, which is one possible quantum analogue of the 
classical class NP ^H] can be obtained by modifying Definition |21 "slightly" . 
Roughly speaking, we change the problem by not asking whether a given 
state |x, 0) is accepted by the circuit Y r in the sense that the output qubit is 



Of C qr (T). 



C qr {T) 




(16) 
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with high probability |1) after applying Y r to the state. Instead, the problem 
is whether there exists a basis state that is accepted. In this case, this basis 
state is called a "witness" for the decision problem. The input string x does 
not define the input state for the circuit. Instead, it defines the circuit Y x 
itself. More precisely, we have: 

Definition 5 (Complexity Class QCMA) 

A language L is in QCMA if there is a uniformly generated family (Y x ) x 
of circuits (where x runs over the problem instances) acting on r + I G 
0(poly(|x|)) qubits such that the following statements are true: 

If we define for every r + l-qubit basis state |y, 0) the r + l — 1-qubit states 
|V>x, y ,o) and |^x, y ,i) by 

Y x \y,0) = a x , y ,o|0) ® |^ x , y ,o) + «x, y ,i|l) ® IVw) , 
then we have: 

1. For every xGl there is a "witness" y such that 

|ax, y ,i| > 2/3. 

2. For every x G" L one has 

|ttx,y,l| < 1/3 

for every y G {0, l} r . 
Then we find: 

Theorem 4 (QCMA-Complete Mixing Problem) 

The following problem is QCMA- complete: 

Given a regular sparse graph G on 2N nodes such that there exists an au- 
tomorphism of G that exchanges all nodes 2j with 2j + l for j — 0,1, . . . , N — 
1. Let ji > and b < a < 1 be two positive numbers with 1/a, 1/6, 1/ (a — b) e 
polylog(iV) and T G polylog(iV) be some time instant. 

Moreover, let N < N be given with the promise that for all j < N we 
have C2j,2j+i(t) G 0(exp(— fit)). Decide if either 

1. there is a number j < N with 

c 2i) 2i + iCr)>ae-" T 

or 



21 



2. for all j < N we have 

C2j,2 j+ i{T) < be-» T . 

To show that the problem is QCMA-hard we construct the graph G x 
corresponding to a given F x according to Section 01 but keep in mind that 
the input string for the circuit (i.e., the potential witness) is here called |y, 0) 
instead of |x, 0). We reorder the nodes such that the states \2j) correspond to 
the states |0)<g>|y,0) and the states \2j + l) to |l)<g> |y, 0) for j = 0, 1, . . . , N, 
where N is chosen appropriately. Section El describes in detail how (1.) and 
(2.) in the above theorem correspond to the cases that the state |y,0) is 
accepted or not by the circuit. Hence we have for each x a graph such that 
a pair (2j, 2j + 1) of nodes exists with j < N and slow decay of probability 
differences if and only if there exists a witness y for the QCMA-problem. 

To see that the problem is in QCMA we recall from the preceding section 
that we can efficiently determine C2j t 2j+i (t) on the appropriate scale by esti- 
mating appropriate expectation values when applying the phase estimation 
procedure to the corresponding states 

-L(\2j)-\2j + l)). (17) 

The numbers j are therefore the potential witnesss: x is in L if and only if 
there exists a j which leads to slow decay of probability differences. We have 
now to construct a circuit where the values j define the possible input basis 
states. 

We define a quantum register with [log A] + 1 qubits where we identify 
each state (fT7|) with 

^(|0>-|1»®L7>- 

Assume we have already constructed a quantum circuit on [logjV] + 1 + 
r qubits (with some appropriate r) describing the whole phase estimation 
procedure including the post-processing of the outcomes with the function 
/ and the repeated sampling. This construction is rather technical, but in 
principle straightforward. Now we modify the circuit such that it accepts 
only those j with j < N. Finally, we add a Hadamard gate acting on the 
left component such that the state can be initialized to |0) instead of the 
superposition |0) — We define the correspondence 

|0)®|i)®|0---0) = |y,0), 
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where j is now allowed to run over all values 0, 1, ... , 2^ logN ~^ instead of only 
0,1, . . . , N . We have now obtained a circuit for which there exists a witness 
y (defining r bits of the r + /-qubit basis state as input for Y x ) iff there exists 
a pair (2j, 2j + 1) with j < N that leads to slow decay. 

7 Conclusions 

We have constructed BQP-complete problems concerning the mixing prop- 
erties of classical random walks. Roughly speaking, the problems are to 
estimate how fast the difference between the probability of being at different 
nodes decays in discrete and in continuous-time random walks. Given that 
the quantum computer is more powerful than the classical computer in the 
sense that BPP ^ BQP we have hence shown that the quantum computer is 
also more powerful in analyzing certain mixing properties of classical random 
walks. 
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